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ABSTRACT 



Context. Mid-infrared spectroscopy of dense illuminated ridges (or photodissociation regions, PDRs) suggests dust evolution. Such 
evolution must be reflected in the gas physical properties through processes like photo-electric heating or H2 formation. 
Aims. With Spitzer Infrared Spectrograph (IRS) and ISOCAM data, we study the mid-IR emission of closeby, well known PDRs. 
Focusing on the band and continuum dust emissions, we follow their relative contributions and analyze their variations in terms of 
abundance of dust populations. 

Methods. In order to disentangle dust evolution and excitation effects, we use a dust emission model that we couple to radiative 
transfer. Our dust model reproduces extinction and emission of the standard interstellar medium that we represent with diffuse high 
galactic latitude clouds called Cirrus. We take the properties of dust in Cirrus as a reference to which we compare the dust emission 
from more excited regions, namely the Horsehead and the reflection nebula NGC 2023 North. 

Results. We show that in both regions, radiative transfer effects cannot account for the observed spectral variations. We interpret these 
variations in term of changes of the relative abundance between polycyclic aromatic hydrocarbons (PAHs, mid-IR band carriers) and 
very small grains (VSGs, mid-IR continuum carriers). 

Conclusions. We conclude that the PAH/VSG abundance ratio is 2.4 times smaller at the peak emission of the Horsehead nebula than 
in the Cirrus case. For NGC2023 North where spectral evolution is observed across the northern PDR, we conclude that this ratio is 
~ 5 times lower in the dense, cold zones of the PDR than in its diffuse illuminated part where dust properties seem to be the same as 
in Cirrus. We conclude that dust in PDRs seems to evolve from "dense" to "diffuse" properties at the small spatial scale of the dense 
illuminated ridge. 

Key words. ISM:individual objects: Horsehead, NGC2023 - (ISM:) dust, extinction - Infrared: ISM 



1. Introduction 

Evolution of interstellar dust was first inferred from measure- 
ments and modelling of the ultraviolet(Uy)-v isible extinction 
(Cardelli et al. 1989; Fitzpatrick & Massal [l986. 1988, Il990t 
[Kim et alJll994l) and from the infrared (IR) co l our variations 
in IRAS and COBE data (.Boulaneer et alJll990H Abergel et alJ 
1 1994 119961; iLaureiis et all Il99lh . The "evolution of interstel- 
lar matter is powered by stars which have a profound influ- 
ence on the chemical and energetic balance of the interstellar 
medium. Along this lifecycle, dust grains undergo efficient pro- 
cessing by shocks and stellar photons which is primarily re- 
flecte d in the dust size distribution (see for instance [jonesll 1 99"7l 
12004 . Such evolution is expected to be particularly significant 
between the dense and diffuse galactic environments. For in- 
stance, the presence of small dust particles in the diffuse galactic 
medium is necessary to account for the near-inf rared (NIR) and 
mid-infrared (mid-IR) o bserved emission (e.g. ISeUgrenl[l984t 
iDraine & AndersonI 119851) while in dense molecular environ- 
men t these small particle s must coagulate with the larger species 
(e.g. lStepnik et alj2003h . Dust has a strong impact on interstellar 
gas because of efficient coupling such as gas heating by the pho- 
toelectric effect or the formation of H2 on grain surfaces. In ad- 
dition, these gas-grain processes are dominated by the contribu- 



tion of small dust grains of radius a < 100 nm (e. g.lHabart et al.l 
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I2OOH IWeingartner & Drainell2001bl: iHabart et al.l |200?). 

Located at the illuminated edges of molecular clouds, pho- 
todissociation regions (here after PDRs; for a review see e.g. 
iHollenbach & Tielenslfl997h are regions of transition between 
dense and diffuse media. In addition to the generally large den- 
sity gradient observed in such regions, there is a gradient of 
the UV-visible radiation field intensity due to the proximity to 
a (some) young star(s) and to the absorption of the radiation 
field by the dust. Consequently, PDRs are the site of signifi- 
cant evolution of both the dust and gas properties which are 
closely interlinked. Dust evolution may be a key to the under- 
stan ding of the unexpectedly large amount of warm H2 in PDRs 
(e.g. lAllers et al.|[2005l; iHabart et al.ll2008l) . PDRs are responsi- 
ble for the reprocessing of most of the UV-visible radiative en- 
ergy output from s tars which is re-emitte d in infrared-millimeter 
wavelengths (e.g. IHollenbach & Tielensi 19971) and consequently 
are important IR emitters of normal star forming galaxies. It is 
thus necessary to understand in detail the emission coming from 
PDRs in order to interpret IR galaxy spectra. 

Infrared Space Observatory (ISO) observations showed a 
systematic diminution of the emission ratio between aromatic in- 
frared bands (AlBs between 3 and 17 fim) emitted by polycyclic 
aromatic hydrocarbons (PAHs) and the mi d-IR continuum emis- 
sion from dense illun iinated ridges ( e. g. Cesarskv et all I2OOOI: 
lAbergel et al.] |2002|) . iRapacioli et all (l2005l) have interpreted 
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Fig. 1. Left panel: The L1630 molecular cloud as seen in the visible (photographic plate obtained using the UK Schmidt Telescope 
and extracted from the Digitized Sky Survey produced at the Space Telescope Institute). The arrow indicates the du^ection of 
the <T Orionis bina i'y. Right panel: NGC2023 (top) and the Horsehead nebula (bottom) as seen by ISOCA M-LW2 at 5-8.5 urn 
jAbergel et alj|20() 2). The contours show areas observed with Spitzer-IRS as part of the "SPECPDR" program jjoblin et alj|2005h . 
The small and large areas delineated by the white line are those observed with the SL and LL modules, respectively. 



this evolution in terms of the chemical evolution of small car- 
bonaceous particle properties by decomposing ISOCAM/CVF 
spectral cubes (5-16 / /m) w ith the Single Value Decomposition 
method. iBerne et al] (l2007l) have extended this study by using 
the Blind Signal Separation method on Spitzer IRS data. A lim- 
itation of these methods is that they do not take into account the 
radiative transfer of stellar photons which may affect the shape 
of the mid-IR spectrum just as dust properties and abundances 
do. The aim of this paper is to study the dust evolution in two 
well known PDRs as seen by Spitzer/IRS and ISOCAM through 
the modelling of the dust mid-IR emission including radiative 
transfer. Although other dust properties like the size distribution 
and/or composition can change, we will interpret observed mid- 
IR spectral variations in term of PAH/VSG relative abundance. 
In section |2] and [3] we present the studied PDRs and the Spitzer 
data. The section |4] is devoted to the observed mid-IR spectral 
shape in these PDRs. In section |5] we present our dust model 
(§ 15.11 ) and the radiative transfer (§ I5.2| | used to model dust emis- 
sion in PDRs. Sections |6] and |2]present results of the modelling 
of the Horsehead Nebula and NGC2023 respectively. We con- 
clude in section[8] 



2. The selected photodissociation regions 

The Horsehead nebula and NGC2023 are located in the LI 630 
molecular cloud as seen in the left panel of Fig. [T] This molecular 



cloud is located at a distance of ~ 400 pc (from the study of the 
distan ces to B stars in the Orion association by Anthony-TwarogI 
119821) in the southern part of the Orion B molecular complex. 

2. 1 . The Horsehead nebula 

In the visible, the Horsehead nebula, also known as B33 
(iBarnardlfTgil) . emerges from the west edge of L1630 as a 
dark cloud hiding the Ha emission of the IC434 HII region (see 
left panel of Fig.[Tli- The Horsehead nebula is a familiar object 
in astronomy and has be en observed many times at visible, IR 
and submm w avelengths dAbergel et alj2003l:lPound et al. 2003 ; 
Tevssier et al.' '2004: 'Petv et al.' '2005^, 'Hilv-Blant et al.l I2OO5I: 
Habartetal. 2005; Ward-Thompson et al. 2006; Phili pp et al. 
2006; Johnstone et al. 2006; Goi coechea et al.ll2006t IPetv et al. 
2007b; Compiegne et al...2007„ for the 21st century only). Both 
the IC434 HII region and the Horsehead nebula PDR are 
mainly excited by cr Orionis which is an 09.5V binary sys- 
tem (Warren & Hesser 1977) with an effective temperature of 
-34 600 K OSchaerer & de Koteiill997 ) and whose direction is 
indicated by the arrow in Fig.[Tl Assuming that cr Orionis and 
the Horsehead are in the same plane perpendicular to the line 
of sight, the distance between them is ~ 3.5 pc (~ 0.5 deg) which 
gives Go ~ 100 (the energy density of the radiation field between 
6 and 13.6 eV in units of the Habing field, Habing 1 968h for the 
radiation field which illuminates the Horsehead nebula. 
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The edge of the Horsehead nebula is delineated by an 
IR e mission filament as seen in the right panel of Fig.[T] 
lAbe rgel et al. (2002) proposed that this filament is due to the si- 
multaneous presence of dense material and a UV radiation field 
from crOrionis. The ISOCAM 5-8.5 yum data seen on the right 
panel of Fig.[T]shows the mid-IR band emission which is known 
to scale with the UV radiation field for Go ranging from 1 to lO'* 
as expected for grains small enoug h (PAHs) to undergo tempera- 
ture spikes dBoulanger et al.lll998l) . Thus, matter located deeper 
in the dense cloud is less excited, which explains this bright nar- 
row IR filament. The Horsehead PDR is an almost plane-parallel 
slab seen edge-on and is thus commonly used as a ben chmark for 
compa rison between models and observations (e.g. iPetv et all 
I2007ah . 

2.2. NGC2023 

NGC2023 is primarily known as a reflection nebula and thus ap- 
pears as an extended source at visible wavelengths (left panel of 
FiglB- The nebula is excited a B 1.5V star (HD37903) which is 
embedded in the L1630 cloud (e.g. Gatlev et al. 1987) and thus 
illuminates the foreground dust (e.g. Burgh et al. 2002) which 
efficiently scatters visible photons (e.g. iDrainei,2003 ). The posi- 
tion of this star, which has an effective temperature of 23700 K, 
is shown on FigIT] Its d istance has been estimated by Hipparcos 
(iPerryman et al.l1l997h to be 470±290pc. The right panel of 
FigHshows NGC2023 as seen by the LW2 channel of ISOCAM 
at 5-8.5 yum which traces mostly AIBs emitted by PAHs. 

The bubble geometry of NGC2023 reflects the influence 
of the young star which blows up the surr ounding matter and 
carves a cavit y in the progenitor cloud (e.g. ICastor et al.ll 19751: 
iFrever et al.l2 003). This scenario i s supported by the l ow extinc- 
tion toward HD37903 (Ay ~ 1.4, 'Burgh etaV lQOl) while the 
total extinction through L1630 is Ay -40 (Hai'vey et al. 1980). 
As seen in FiglT] NGC2023 is thus basically a spherical cavity 
seen face-on. 

In this work, we focus on NGC2023 North (hereafter 
NGC2023N) which exhibits a strong variation of its mid-IR 
spectrum as shown bv .Abergel et al., (.2002) . 



3. Observations 

3.1. Infrared spectrograph data 

The Horsehead Nebula and NGC2023N were observed as part of 
the "SPECPDR" program (Jobhn et aj1l2005l) . using the infrai-ed 
spectrograph (IRS, Houck et al. 200 3) whi ch is a slit spectro- 
graph on board Spitzer (We rner et al 120043) . The present obser- 
vations were performed with the spectral mapping mode and us- 
ing the short-low (SL, 5.2-14.5yum, slit size: 57" x 3.6", R=64- 
128) and long-low (LL, 14-38 yum, slit size: 168" x 10.6", R=64- 
128) modules of the instrument. The integration times were 14 
and 60s per pointing for the second (5.2-8.7/im) and the first 
(7.4-14.5;um) orders of SL, respectively, and 14 s per pointing 
for both orders of LL. 

Starting from BCD products, we have developed a pipeline 
which builds spectral cubes (two spatial dimensions and one 
spectral dimension) in a homogeneous way from the data de- 
livered by the S pitzer Science Center (SS C). This pipeline has 
been presented in lCompiegne et aLl(l2007l) . The LL data between 
35 and 38 fim are not considered due to the strong decrease of 
sensitivity at these wavelengths. We therefore obtain 5 to 35 fim 
spectral cubes for the observed area. The absolute photometric 
uncertainty is ~10%. 
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Fig. 2. Mean IRS-SL and IRS-LL spectra of NGC2023N and 
the Horsehead Nebula. The dot-dashed lines indicate the wave- 
lengths of the main H2 rotational lines 0-0 S(4) to S(0) at 
6.9, 9. 7, 12.3, 17.0 and 28.2 //m which are analyzed else- 
where (Habar t et 31.11200 8*). The dotted lines on the Horsehead 
spectrum indicate the wavelengths of the fine structure lines 
[ArII]6.98;um, [Nell] 12. 8 //m, [SIII]18.7 and 33.4;um, [Sill] 
34.8yum. 



The zodiacal emission Z(A) must be subtracted from these 
data. We assume that this contribution does not vary across 
the observed area of the sky for a given object. We estimate 
Z(A) using the SSC backgro und estimator ] which is based 

etan iT 



on the COBE/D IRBE model dKelsaU 
^Com piegne et alj (|2007) , we take into 



Il998h . Following 
account the zodiacal 
contribution in the dark observation (at position RA - 268 ^96, 
DEC = 65 °43 in the sky) in the BCD level of the IRS data. 



3.2. ISOCAM data 

Both the Horsehead nebula and NGC2023 were observed with 
ISOCAM on board the Infrared Space Observatory (ISO) in the 
LW2 (5-8.5 jum) and LW3 (12-18//m) broad band photometers 
(lAbergel et al.ll2002l) . The data shown in Fig. [I | and Fig. [3] has 
been reduced following lMiville-Deschenes et al.l (12000). For the 
two filters, the photometric accuracy for the extended emission 
is ~ 10%. The zodiacal contri bution subtrac ted from this data has 
also been estimated from the lKelsall et al.l (I1998i) model. 



See |http://ssc.spitzer.caltech.edu/documents/background/| 
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Fig. 3. Map of the LW2/LW3 ratio in NGC2023. The contours are the H2 v= 1-0 S(l) Hne at 2.12;um data as seen with the SOFI 
instrument at the ESO New Technology Telescope. The spatial resolution of SOFI data and of LW2 were degraded (with a Gaussian 
convolution) to match that of LW3 (~4.5")- The star indicates the position of HD37903. The two circles indicate the locations of the 
spectra seen in the right panel. The upper spectrum is from the northern position (deep in the cloud) and the lower spectrum is from 
the southern position (close to the star). The 5-14.5 jjm part of the spectra is IRS-SL for the upper spectrum and ISOCAM-CVF for 
the lower spectrum (and IRS-LL at longer wavelengths in both cases). 



NGC2023N has also been observed with ISOCAM-CVF 
dAberpelet al.ll2002l) . The area observed with ISOCAM-CVF 
includes almost entirely that observed with IRS-LL and com- 
pletely that observed with IRS-SL. We use this data to obtain the 
same spatial coverage for spectral mapping data for the 5-16/im 
range as for the 14-35yL(m range (IRS-LL). The spectral resolu- 
tion of ISOCAM-CVF is A/AA -40. The data reduction and 
zodiacal light subtraction is described in lBoulanger et al.l (l2005h 
and the corresponding photometric accuracy is ~ 20%. 

4. Spectral variations 

Mean mid-IR spectra of both objects are displayed in Fig. [31 
They are typical of PDRs and present emission in the AIBs 
at 6.2, 7.7, 8.6, 11.3, 12.7 fj.m, the mid-IR dust continuum as 
well as in H2 rotational lines within the first vibrational level. 
The NGC2023N spectrum, which is brighter, also clearly ex- 
hibits the 16.4yL(m band and the 17.1 and 19 u m bands already 
reported in NGC7023 bv lWerner etsl\ (l2004bl) . The spectrum of 
the Horsehead nebula also presents fine structure lines of ionised 
species (dotted lines in Fig.|2]l which arise from the HII region 
surrounding the Horsehead nebula. In this pa per, we hmit our 
study to the neutral gas zone (i.e. the PDR, see iCompiegne et al] 
120071 for study of the HII region). 

Fig. [3] shows the map of the band ratio LW2(5- 
8.5 um)/LW3(12- 1 8jUm) in NGC2023. As akeady noted 
by lAbergel et al.l (|2002|) . this ratio varies across dense, pho- 
todissociated interfaces. The northern PDR is well delineated 
by the H2 v- 1-0 S(l) emission along a nearly east-west ridge 
at ^2000 ~2°13' (between the two circle positions in Fig.[3j 
which spatially correspond to a striking variation of LW2/LW3. 
This ratio indeed decreases by a factor of ~2.6 between the 
two positions marked by circles in Fig.|3] for which spectra 
are displayed in the same figure. From their analysis of the 
ISOCAM-CVF (5-16//m) spectral maps, lAbergel et al.1 (|2002|) 
concluded that this evolution of LW2/LW3 is related to a 
decrease of the AIBs relative to the mid-IR continuum emission 
while moving away from the star and into the molecular cloud. 



This conclusion is clearly confirmed by the present IRS spectra 
which also extend to longer wavelengths and better define the 
mid-IR dust emission. 

The spectra of Fig.|3]also show a variation of the shape of the 
20-35 jum continuum between the two positions. While almost 
linear deep into the cloud, it shows a non-linear rise closer to the 
star Note that bigger grains which are at thermal equilibrium 
could be hot enough to emit below /I ^ 35 yum for positions close 
to the star while they are colder deeper in the cloud. 

Conversely, the Horsehead observations do not show signif- 
icant spectral variations. This may be due to the sharpness of 
the photodissociated interface barely resolved by IR data and 
across which the variations occur However, we note that while 
in NGC2023N the LW2/LW3 ratio goes from -1.7 close to 
the star to ~0.6 deeper into the cloud, the value of this ratio 
for the unresolved filament of the Horsehead nebula is ~1.1. 
This value is almost median between the extreme values across 
NGC2023N. In fact, the mean spectra of Fig.|2] have similar 
LW2/LW3 values (~ 1) suggesting that we see "average" prop- 
erties in the Horsehead. 

5. Dust emission model for PDRs 

In order to analyze the IR spectra of PDRs and to trace dust 
abundance variations by taking into account excitation effects, 
we developed a dust emission model and also treat the radiative 
transfer problem inside a PDR. Our dust model is an updated 
version of that Desert et al. (1990) (hereafter DBP). We combine 
this model with a 1 dimension radiation transfer model which is 
fully consistent with the dust properties of the dust emissivity 
model. Modifications implemented to the DBP model are de- 
scribed in section lSTI while the description of the transfer model 
is given in section lS^ 

5. 1 . The dust model 

The interpretation of dust observables (extinction and IR emis- 
sion) in the ISM requires at least three dust populations of in- 
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Table 1. Parameters of the AIBs with Nc and Nh the num- 
ber of carbon and hydrogen atoms in the molecule respectively. 
CC(CH) bands strength scale with Nc (Nh)- 



creasing sizes (e.g. DBP. lLi & Drainell200ll:IZubko et al.ll2004l) . 
In DBF, these are: (i) large (radius a ~ 0.4 to 1 nm) polycyclic 
aromatic hydrocarbons (PAHs) which account for the AIBs and 
for the far-ultraviolet (FUV) non-linear rise in the extinction 
curve, (ii) carbonaceous nanoparticles (a ~ 1 to 10 nm) called 
very small grains (VSGs) which carry the mid-IR continuum 
emission and the extinction bump at 2175 A and (iii) big grains 
(BGs, a ~ 10 to 100 nm) which produce the far IR emission, 
the l/A rise at visible near-IR wavelengths as well as the 10 and 
20 fim absorption bands in the extinction. In the diffuse interstel- 
lar medium, small grains (PAHs and VSGs) undergo temperature 
spikes triggered by the absorption of stellar photons and cool by 
emission in the near and mid-IR range. Conversely, BGs, which 
have a longer cooling time and a shorter timescale between ab- 
sorption of two photons because of their size, stay at constant 
temperature and emit like grey bodies. 

The dust properties that we describe below allow the re- 
production of the extinction curve and emission spectrum of 
the diffuse high galactic latitude ISM (so called "Cirrus") as 
seen in Fig.|4] The exciting radiation is the Interstellar Standard 
Radiation Field (ISRF, Mathis et al. 1983). We call this dust pop- 
ulation Cirrus dust as a reference to compare to dust emission in 
other interstellar environments. 

The absorption cross-section of PAHs has been updated. 
We tak e the visible-UV cross-section from 'V erstraete & Legerl 
(1 1992b and apply their size dependent cut-off for electronic tran- 
sitions in the near-IR. The resulting cross-section provides a 
good match to the available laborator y measuremen ts on species 
containing 20 to 30 C atoms (Jobl in et al.l |I992|) . The cross- 
sections of the AIBs (a t 3.3, 6 . 2, 7.7 8.6, 1 1.3, 12.7 and 16.4 //m) 
are based on Verstraete etal] (|200l)) and iP ech et"aL l (|2002|) and 
assumes singly ionized PAHs. ISO-SWS spectra also indicates 
the presence of a broad band at 6.9 ^m and a band at 12yum 
(I Verstraete et al.ll200ll; iPeeters et aT] |2002). In addition, Spitzer 
data has shown the presence of another band at 17 /im w hich be- 
haves like the rest of the AIBs (Peeters et ak 2004; Smith et al.l 
I2007h . We therefore include these bands in our AIB list which is 
summarized in Table[T] The positions and widths of these bands 
were derived from ISO and Spitzer spectra. Peak values of cross 
sections of all AIBs are adjusted in order to match the Cirrus 
observed spectrum. We assume planar and fully hydrogenated 
PAHs with hexagonal symmetry (ZJgft). The molecular radius a 
is related to the number of carbon atoms per molecule Nc as 
a(A) = 0.9 VA^c the number of hydrogen atoms is given 
by Nh = "s/^ Nc- Finally we take the far-IR cross section from 
ISchutte et al.l (119931) . 
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Table 2. Parameters of the size distribution n(a) for the three 



components of our dust model. We assume that n(fl) 
a between amin and oimax- 



with 



VSGs are assumed to be graphitic and both the optical prop- 
erties and heat capacity have been updated (e.g. I Praine & Lil 
1200 1 1) . Their optical properties have been calculated with Mie 
theorjQ applied to small graphitic spheres in the 2/3- 1/3 ap- 
proximation (e. g. .Draine & Lee 198 4). The dielectric constant 
was taken from lLaor & Drainel(ll993b . 

As in DBP, BGs are taken to be silicates with carbonaceous 
mantles or inclusions in order to provide the required absorption 
efficiency in the NIR compared to bare silicates. The near-IR ab- 
sorption cross section of BGs now has a /I ' ^ dependence which 
allows a better match of the observed extinction curve. 

We assumed the size distributions of all the dust popula- 
tions of our model to be a power-law type, n{d) oc a'^-^ with 
minimum and maximum radii flmi,, and fl^ax listed in Table|2] 
This simple form of size distribution original ly introduced by 



Mathis et al. ( 197 7J) re r nains a good paradigrn ( Kim et aLlll994i 
Weingart ner & Drain3 l2001al: IClavton et al.1 12003). Note that 
abundances of the three populations remain unchanged from 
those of the DBP model. 

One of the consequences of the update is that UV-visible 
cross sections (see Fig.|4]i of the PAHs and VSGs are more sim- 
ilar than in the DBP model. In fact, these two species are as- 
sumed to contain mostly sp^ hybridization carbons which are 
responsible for the 2175 A bump (n — > n* electronic transition) 
and for the non-linear rise in the far UV (bump at 800 A due 
to a o" ^ cr* electronic transition). This similarity has conse- 
quences for the possible self shielding of the two species which 
will cause an evolution of the spectral shape depending the shape 
of the exciting radiation field spectrum. 



5.2. Radiative transfer 

We describe here how the total exciting flux F'°'(z, /I) is esti- 
mated at each depth z of the cloud. This total exciting flux results 
from the contribution of the illuminating radiation field trans- 
fered into the cloud Fiz, A) and of the dust emission Fdust(z, A), 
as illustrated by Fig|5] 

We first transfer the incident radiation field into the cloud 
which is represented by a semi-finite plane-parallel slab of pro- 
ton density nniz) at each depth z (z = at the free interface). We 
use the dust properties described above. Note that the PAHs are 
assumed to have a zero albedo and that the albedo of VSGs is 
low due to their small size. The scattering cross section is thus 
dominated by the contribution of BGs. 

The incident radiation field Fq (X) is the sum of the mean in- 
terstellar radiation field (ISRF) of iMafliis et al.1 (Il983h and of a 
blackbody whose flux is diluted by {Ri^jdf' with d the cuiTent 
distance to the star and 7?* the stellar radius. Note that the flux of 
the ISRF is obtained by integrating the intensity over 2n stera- 
dians since the cloud is illuminated on one side. Moreover, we 



Code retrieved from 'h ttp://atol.ucsd.edu/scatlib/index.html 
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Fig. 4. Upper panel: left and right: Diffuse interstellar medium extinction. The solid line shows our model. Diamonds are the 
mean extinction curve from Savage & Mathis ( 1979) while triangles are the mean extinction curve of Mathis ( 1990). Lower panel: 
left Emission spectrum for the diffuse interstellar medium for Nh - 10^*^ cn T^. The sohd Une is our dust model. Data points are 
(i) the 3.3 jum emission (triangle) of cirrus associated with the molecular ring ("Giard et aLlll994}) . (ii) the ISOCAM-CVF spectrum 
between 5 and 16^m of the diffuse Galactic emission at b~l°(Flagev et al. 2006), (i ii) the DIRBE measurements (squares) for 
\b\ > 25° (lArendt et alj|l998l) and (iv) the FIRAS submm specti'um of Boulanger et ail (1199 6V Right: Diffuse ISM model emission 
spectrum (thin solid line) over the wavelength range of IRS. The ISOCAM-CVF spectrum is superimposed (thick solid hne). The 
model spectrum has been smoothed to the resolving power (R~40) of ISOCAM for A < \5 //m. 



limit the radiation field to photon energies of less than 13.6 eV 
since we consider neutral gas (PDR). 

The plane-parallel approximation allows us to simplify the 
treatment of the scattering. We assume that each incident pho- 
ton is either scattered forward deeper into the cloud (increasing 
z) or scattered back towards the cloud interface (decreasing z). 
Assuming that a photon can only be backscattered once, the flux 
of the radiation field at a given depth z into the cloud can be 
written as 



F{z,A) = F,(z,A) + Fh{z,Zm...X) 



(1) 



where F,(z, A) is the transmitted flux from the illuminated sur- 
face of the cloud and Fh{z, Zmax, ^i) the flux backscattered by the 
part of the cloud deeper than z- We impose a finite depth Zmax for 
the cloud in order to estimate the backscattered flux. 

With our assumptions, we only need absorption cross sec- 
tions (cr„(/l)) and backscattering cross sections (cri,(A)) to per- 
form the transfer since photons that are scattered forward are 
simply treated as transmitted. Our dust model only computes the 
integrated scattering {cTsiA)) cross section and not the differential 



one {dcTsiA) / dO). We thus use the lHenvev & Green steinid 19411) 
phase function (Oh-g(6')) to calculate the backscattering cross 
section as 



(TbiA) - 2jio-s(A) 



Jnll 



sinO do. 



(2) 



Note that when taking g = < cosiff) > - 0.6 which is the 
common value used, in ac cordance with the diffuse interstellar 
medium observations (e.g. I Witt et alJ[T997l : ISchiminovich et al.l 
1200 ll) . ~ 12 % of scattered flux is backscattered. 
The transmitted flux is given by 

F,(z,A) = Fo(^)e-"»'("^-''> (3) 

with the opacity Text(0 — > z. A) defined as 

Text(zi ^Z2,^) = £ nH(z')(cra(z',A) + crb(z',A))dz' (4) 
The backscattered flux Fi,iz, Zmax, ^) is given by 



Fbiz,Zm!ix,yl) = dFbiz,A)e 



-T,„(Z'^Z,-1) 



(5) 
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Fig. 5. Schematic view of the plane-parallel geometry used. 

A) is the total exciting flux which is considered in deriv- 
ing the dust emissivity. Fdust(z, ^) is the flux related to the dust 
emission from the cloud itself. Ft{z, A) and Fi,{z, Zmax, ^i) are the 
flux related to the illuminating radiation field transmitted from 
the illuminated edge and backscattered by the deeper part of the 
cloud, respectively. 



with 



dFb(z',A) ^ F,(z', A)dTb(z\ A) 



(6) 



where dTh(z', A) is the backscattering optical depth of the layer 
dz' at a depth z'- 

From the exciting flux F(z, A) due the illuminating radiation 
field, we derive the dust emissivity e'{z. A) which will be used to 
calculate the exciting flux /^du,st(z, ^) due to the dust emission. In 
the plane-parallel approximation, we can estimate this contribu- 
tion as 



Fdustiz,A) = i 



nniz') e'iz. A) e 



dz' 



(7) 



with e'(z'. A) the dust emissivity resulting from the illuminating 
radiation field F(z, A), nniz') the density profile and Text(z' — » 
z. A) the extinction opacity between the depth z' and z as defined 
by eq.|4] The 1/2 factor is related to the fact that the dust emis- 
sion is assumed to be isotropic and that consequently, half the 
power emitted at a depth z' does not reach the depth z, since the 
dust scattering cross section is very low at the relevant wave- 
lengths. Note that for the studied depth (Ay :$ 20), an iterative 
computation of the dust emission contribution to the dust heating 
is not required since the mid and near IR dust emission, which 
contributes significantly to the dust heating is mostly emitted by 
the dust at low depth, whose heating is completely dominated by 
the UV-visible photons of the illuminating radiation field. 

We can finally derive the dust emissivity e{z, A) as a function 
of the depth in the cloud, as excited by A) - F(z, A) + 

FiMstiz, A). 



6. Modelling the Horsehead nebula 

Due to its almost edge-on geometry, the Horsehead nebula can 
be modelled using plane-parallel geometry. 

The radiation field is characterized by - 8.5/?© 

and Teff — 34620 A" corresponding to an 09.5V star 
dSchaerer & de Koteii[T997h located at a distance o f 3.5 pc from 
the cloud. We use the density profile deduced by iHabart et alj 



(|2005|)from a comparison of model gas emission profiles to ob- 
servations of H2, CO and dust emission: 



(8) 



with 



210^ 



.-3 



0.02 pc (10") andy6=4. We inte- 



grate the output emissivity e{A, z) along the line of sight to obtain 
intensities in units of MJy sr"', those of our data. We further im- 
pose a length for the PDR along the line of sight, namely, /pdr. 
We neglect the absorption at these wavelengths since the visual 
extinction range between Ay ~ 2 - 30 as function of z in the 
line of sight along the Horsehead. As in Habart et alj (I2005h . we 
perform a 6° rotation of the PDR to the line of sight. 

6. 1 . Emission profiles 

Fig. |6] shows observed and modelled profiles at 7-9 //m and 22- 
24yum normalized to the peak value. These spectral ranges are 
unambiguously related to the AIBs (7-9 yum) and the mid-IR con- 
tinuum (22-24 yum) which was not the case for the ISOCAM- 
LW3 (12-18//m) broad band which contains both mid-IR contin- 
uum and certain AIBs (the 12.7 16.4, 17.1 fim). Modelled pro- 
files were convolved by a Gaussian whose width corresponds 
to the IRS spatial resolution at the corresponding wavelengths 
(1.9"at 1-9 jjm and 5.5"at 22-24jum). The observed profile is 
defined as the mean over the strip shown in Fig.|6] Both the ob- 
served and modelled angular offsets are defined with respect to 
the peak emission position at 7-9 fj.m. We can then see that the 
mid-IR continuum (22-24 jum) and AIB (7-9 fim) emission peak 
at the same depth in the cloud and are well reproduced by the 
model. 

Although we take into account the profile broadening caused 
by the point spread function and by the tilt of the illumi- 
nated ridge, we see that the modelled profiles are narrower than 
the observed one. This must be related to differences between 
our plane-parallel representation and the real geometry of the 
Horsehead nebula (e.g. multiple filaments superimposed). The 
absolute intensity of the mid-IR continuum (22-24 jum) at the 
peak position can be reproduced if /pdr ^O.OSpc. Th i s valu e 
is only 20% lower than that found by iHabart et alj (l2005h . 
However, the modelled AIB intensity is then 2.4 times higher 
than the observed one (i.e would have required /pdr = 0.033 pc 
for the absolute intensity to be reproduced). We can then con- 
clude that taking the radiative transfer effect into account, 
"Cirrus dust" properties cannot explain the observed AIBs/mid- 
IR continuum ratio. 



6.2. Spectral comparison 

The AIBs excess with respect to the mid-IR continuum can also 
be seen in Fig.|2] In both model and observations, spectra of this 
figure are defined as the average over positions whose intensities 
at 7-9 fim are at least 2/3 of the peak values. For the observed 
spectra, we must ensure that we are looking at the same position 
within the cloud (i.e. same depth) since we are looking at wave- 
lengths ranging from 5 to 33 fim whose spatial resolution are not 
the same. We therefore bring the resolution at all wavelengths 
to 7.9", i.e., that of the 33 fim by convolving the spectra with 
a Gaussian of appropriate width. For the modelled spectra, we 
simply smooth with a Gaussian of 7.9"width. As already seen 
above with the 1-9 fim and 22-24 //m profile comparison, the 
modelled spectrum shows an excess of AIB emission relatively 
to the mid-IR continuum. Furthermore, the model shows that the 
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Fig. 6. Left panel: Image of the Horsehead Nebula has seen by IRS-LL. The two cuts show the limit between which the mean 
observed profiles are computed. Right panel: Comparison of the observed and modelled profiles at 7-9 //m (AIBs) and 22-24 yum 
(mid-IR continuum) normalized to the peak values. The offset origin of the cuts correspond to the 7-9 fim emission peak. crOrionis 
is located toward negative offsets. 




Fig. 7. Comparison between observed (thick line) and modelled (thin line) spectra at the peak position of the emission of the 
Horsehead Nebula (see § l6.2l for the exact definition). In the left panel, we show Cirrus dust emssion. In the right panel, the PAH 
abundance has been reduced by 2.4 with respect to the Cirrus abundance. Note that IRS spectra have been cut at A ~ 33 //m to 
avoid the fine-structure lines [SIII] 33.4 jum and [Sill] 34.8 fim. 



mid-IR continuum in the IRS range is dominated by VSG emis- 
sion. We can therefore discuss the AIBs to mid-IR continuum 
ratio in terms of PAH/VSG relative abundance. 

Fig. |7] also shows the case where the PAH abundance is re- 
duced by a factor of 2.4 with respect to that of CiiTus: the match 
to the data then becomes very good. We thus conclude that the 
PAH abundance relative to VSGs is 2.4 times higher in the CiiTus 
medium than at the peak emission position of the Horsehead 
Nebula. 

6.3. Effect of differential excitation 

Here we discuss and quantify the effect on the mid-IR emission 
spectrum of differential excitation of PAHs and VSGs due to ra- 
diative transfer. As discussed in § 15.11 differences in the absorp- 
tion cross sections can have an impact on the relative excitation 
of different dust populations depending on the spectral shape of 
the exciting radiation field. 

Fig.[8]shows the modelled 7-9 /^m/ 22-24 //m (i.e AIBs / mid- 
IR continuum) ratio evolution as the function of the optical 



3.0 r 1 1 1 — I 1 1 1 — r|3.0 

I„(7-9/.im)/I„(2a-24/im) - 

Pass, pah / ^abs, vsc 




0.0 [ I I I I I I I I I I I \ I III lI 0.5 

0.1 1.0 10.0 

Av 

Fig. 8. The 7-9 /vm / 22-24 //m emission ratio and the ratio of 
absorbed power per C atom by PAHs and VSGs as a function of 
the optical depth into the Horsehead Nebula. 



M. Compiegne et al.: Dust processing in photodissociation regions 



9 



1 — I — I — I I 1 1 1 




10"' I I I I I I I I I 

0.1 1.0 10.0 100.0 

Fig. 9. Upper panel: Exciting radiation field F^°\z,A) for 
Ay = 0.1, 1, 4, 7, 10, 20 with increasing thickness of the line for 
increasing Ay and (lower panel) the absorption cross sections 
per C atom of PAHs and VSGs of the model (1 Mb = 10 cm^). 



depth into the cloud for the Horsehead Nebula model as describe 
above. The dust properties used are those of Cirrus. The col- 
umn density (outpu t of the model) was converted to Ay with the 
iBohlin et alj 11978 1) ratio, Nn/Ay ~ 1 .87 lO^' cm-^ mag-\ We 
see that the ratio increases by ~ 25% at 0.1 ^ Ay ^ 3 and then 
decreases by a factor of ~4.5 at Ay ^ 10. Fig.[8]also shows that 
the 7-9 jum / 22-24 /^m ratio evolution is related to the evolution 
of the absorbed power by PAHs and VSGs (per C atom) as a 
function of depth. Note that differences between the two curves 
are related to the fact that variations of the temperature distribu- 
tion of PAHs and VSGs (which involve stochastic heating) as a 
function of the optical depth also cause spectral shape evolution 
of the individual population emission spectra. 

Qualitatively, this evolution of the absorbed power ratio can 
be understood by looking at the spectral evolution of the excit- 
ing radiation field F^°\z, A) (upper panel Fig.|9]l and the cross 
sections of PAHs and VSGs per C atom (lower panel Fig.|9]). 
The absorption cross section of PAHs is greater than those of 
VSGs for A :$ 1 .5 yum and lower for 1.5 ^ A i 6 fim. Hence, 
the reddening of the exciting radiation field caused by extinction 
successively favours PAH and VSG excitation. 

We emphasize that the observed AIBs to mid-IR continuum 
ratio of the Horsehead nebula cannot be explained by such dif- 
ferential excitation. In fact, the 7-9 jum / 22-24 yum ratio is equal 
to 0.9 for the peak emission spectra (Fig.|2ll, compared to the 



Cirrus value of 2.1. In the model, the peak emission occurs at 
0.2 :$ Ay :$ 0.8 (where the profile at 7-9 //m is equal to 2/3 
of the peak value following the definition of the peak spectrum). 
For this depth, the ratio does not vary by more than 13% of the 
CiiTus value. To account for the observed variations with only 
the radiative transfer effect, the mid-IR spectrum at the peak po- 
sition should be emitted by dust at a depth of Ay ~7, corre- 
sponding to a decrease of more than three orders of magnitude 
of the dust mid-IR emissivity compared to Ay =0. 

7. Modelling of NGC2023N 

As already pointed out in section|4](see also lAbergel et alj2002l) . 
unlike for the Horsehead nebula, we spatially resolve the spec- 
tral shape evolution in NGC2023N. For the two spectra of Fig. [3] 
the 7-9 yum / 22-24 fj.m ratio goes from 1 .9 close to the star (lower 
spectrum of Fig. [3]) to 0.4 for the northern position (upper spec- 
trum of Fig. [3]). The former must be emitted by a relatively dif- 
fuse medium around the star just at the external edge of the dense 
illuminated ridge traced by the H2 v= 1-0 S(l) emission while 
the latter is emitted by denser gas at the ridge. 

7. 1. The cavity spectrum 

It is interesting to note that the 7-9 /vm / 22-24 yum ratio of the 
spectrum from the diffuse illuminated part of the PDR (~ 1.9) is 
close to the Cirrus value (~2.1). As shown in Fig.[TOl the ob- 
served spectra can be reproduced by the model by using Cirrus 
dust properties. The modelled spectrum was obtained by con- 
sidering a star of T,ff = 23700 K and R^^eRe (for B1.5V star) 
located at 0.3 pc which corresponds to the ~150" between the 
star and the location of the observed spectrum for a distance 
of 400 pc. The slope of the continuum emission is reproduced 
by applying an extinction of Ay =1.25 between the star and 
the emitting dust. Less extinction would result in hotter BGs 
and in a too steep slope of the continuum while more extinc- 
tion would result in colder BGs and a slope not steep enough. 
That extinction is consistent with the presence of matter in the 
cavity (|Witt et al. 1984; Burgh et al. 2002). Moreover, this value 
of Ay ~ 1.25 caused by the cavity matter is of the same order 
as tha t reported between the Earth and the star by iBurgh et"al] 
(I2002h . which is Ay ~ 1.4. It is not solely the extinction that al- 
lows us to adjust the BG emission. Considering the error bars on 
the distance estimate of NGC2023, an underestimation of this 
distance would cause an underestimation of the applied dilution 
of the star radiation field in the model. The absolute intensity of 
the modelled spectrum matches that of the observed one for a 
column density of Nh = 4.6 10^' cm^^ (in the line of sight). 

As suggested in section|4] we see that the non-linear rise of 
the mid-IR continuum can be explained by a significant emis- 
sion of BGs at A ^ 25 fim while according to the model the 
continuum emitted by the VSGs is mostly linear. Note that the 
intensity of the exciting radiatio n field after dilution and extinc- 
tion isx ~ 120 (in unit of the Habingj (Il968i) field between 6 and 
13.6eV) which is enough for the BGs to emit noticeably in the 
IRS spectral range. 

BGs do not emit significantly in the 22-24 fj.m range at this 
position and will be colder further from the star which make 
them emit at longer wavelengths. The decrease by a factor of 
~ 5 of the 7-9 jum / 22-24 /urn ratio may then be related to the 
evolution of the PAH/VSG relative emission, as in the Horsehead 
case. This is coherent with the linear shape of the mid-IR contin- 
uum for the observed spectrum deep into the dense ridge (Upper 
spectrum of Fig.[3]l 
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Fig. 10. Comparison between the modelled and observed spectra 
of the external diffuse part of the PDR in NGC2023N. The po- 
sition of the observed spectrum is the same as the lower spectra 
of Fig.[3] The observed spectrum is ISOCAM-CVF at 5-14.5 yum 
and IRS-LL at longer wavelength. The modelled spectrum was 
brought to the ISOCAM-CVF spectral resolution for /I < 15 yum 
and IRS Low-Res for /I > 15 ;um for comparison. The dust prop- 
erties are those of Cirrus. 



AIB spectrum is also well reproduced by the Cirrus RA^H 
properties which suggest similar properties of these emitters in 
the diffuse medium around the star and in Cirrus. 

7.2. The spectral shape evolution 

As seen in Fig|3] the spectral shape evolves in the dense ridge 
traced by the H2 v- 1-0 S(l) emission located just after the lo- 
cation of the studied cavity spectrum. As in the case of the 
Horsehead, we can model the evolution of the 7-9 fim / 22-24 /im 
emission ratio as a function of the optical depth into the dense 
illuminated ridge in order to quantify the radiative transfer ef- 
fects. Unlike the Horsehead nebula, t he density profile i s not 
constrained. However, as reported by iField et aLl (Il998l) (see 
references therein), the density must range between lO'* and 
10^ cm^^ in the dense illuminated ridge as shown through the 
use of steady-state PDR models to study the brighter part of 
NGC2023. The simplest model for the dense ridge is then a den- 
sity wall illuminated by the B1.5V star located at 0.3 pc. The 
radiation field of the star that illuminates the ridge must be ex- 
tinguished with Ay = 1.25 to be consistent with the value of the 
extinction at the location of the studied cavity spectrum, which 
is the position where the dense ridge arises. Note that this ex- 
tinction corresponds to a constant density of «// - 2550 cm^^ 
for the cavity. 

FigfTTI shows the result obtained for the 7-9 /vm / 22-24 yum 
emission ratio as a function of the optical depth in the dense 
ridge for densities of 10"* and 10^ cm~^. We can see that results 
for the two different densities are similar and also are similar to 
the result obtained in the Horsehead case. Hence, radiative trans- 
fer effects cannot explain the spectral variation in NGC2023N 
since a depth of Ay ^ 12 is required to explain the value of 0.4 
for the 7-9 ;um / 22-24 /zm ratio observed in the deep/dense part 
of the ridge. Such an extinction would cause a decrease of more 
than 3 orders of magnitude of the the 22-24 jjm dust emissivity 
compared to the edge of the ridge (Ay - 1.25) while the ob- 
served intensity increases from ~70MJysr"' to ~100MJysr"' 
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Fig. 11. 7-9 yum / 22-24 /vm emission ratio as a function of the 
depth in the dense illuminated ridge of NGC2023N. 



between the cavity spectrum and the deep/dense ridge spec- 
trum. We can conclude that the PAH/VSG relative abundance 
increases by up to a factor of ~ 5 from the deep/dense part to 
the diffuse/illuminated part of NGC2023N. Moreover, dust prop- 
erties in the diffuse medium surrounding the star in NGC2023 
must be the same as in the diffuse high galactic latitude medium 
(Cirrus). 

8. Conclusion 

We presented mid-IR spectral imaging observations of the 
Horsehead nebula and NGC2023N obtained with the infrared 
spectrograph on board the Spitzer Space Telescope. These ob- 
servations allow us to confirm the AlBs / mid-lR contin- 
uum evolution at dense illuminated ridg es already observed 
with the Infrared Sp ace Observatory (e.g. lAbergel et al.l 120021 : 
iRapacioli et"ani2005h . 

We developed a new dust emission model based on the 
iDesert et alj d 19901) model. This model successfully reproduces 
the emission and the extinction curve of the diffuse interstellar 
medium at high galactic latitude (i.e. the Cirrus). The Cirrus dust 
properties are then used as a reference to be compared to dust 
properties in the two studied PDRs. This dust emissivity model 
is coupled to a radiative transfer model in order to properly take 
into account excitation effects on the mid-lR spectral shape evo- 
lution. 

For the PDRs studied here, we conclude that excitation ef- 
fects cannot account for the observed AlBs / mid-lR continuum 
evolution. We interpret this evolution in term of PAH/VSG rel- 
ative abundance variation since these two species dominate the 
5-35 jum spectrum of the two studied PDRs except at A ^25 fim 
for positions close to the star (in the cavity) in NGC2023 where 
the BG continuum is significant. 

In the Horsehead Nebula case, we do not spatially resolve 
any spectral variation. We conclude from the study of 7-9 //m 
(AIBs) and 22-24 jim (mid-IR continuum) emission profiles and 
of the spectrum at the emission peak position that the PAH/VSG 
relative abundance is 2.4 times lower than in the Cirrus. In 
NGC2023N, we resolve the spectral variation. We modelled the 
spectrum of the diffuse medium around the star for a position 
at the external edge of the dense illuminated ridge traced by H2 
v=l-OS(l). This spectrum can be reproduced by using Cirrus 
dust properties. We attributed the decrease of a factor of ~5 of 
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the 7-9 //m/ 22-24 jum (AIBs / mid-IR continuum) from the dif- 
fuse illuminated part to the deep/dense part of the PDR to a de- 
crease of the PAH/VSG relative abundance to ~ 1/5 of the Cirrus 
value. It then seems that dust properties evolve fully "from dense 
to diffuse" properties within the small spatial scale of the dense 
illuminated ridge. Note that the obtained PAH/VSG relative 
abundance for the "unresolved" Horsehead is almost the median 
between the diffuse and dense part of NGC2023N, indicating 
the mean properties of the observed very small particles in the 
Horsehead. 

An increase of the PAH/VSG relative abundance from the 
dense/deep to the diffuse illuminated part of the PDR is fully 
consistent with a sc enario of photoevaporation PAH clusters (i.e , 
of the VSGs) (e^ . ICesarskv et all 120001: iRapacioli etal] 120051; 
iBerne et al]|2007l) . Other physical processes could explain this 
variation, such as (i) the scenario evoked in Jones (2005) of the 
evolution of aliphatic hydrocarbons (which do not emit AIBs) 
which are present in the dens e media (e.g. I Jones et aTl Il990t 
|Jonesl[T990t iDartois et al.ll2007l) into aromat ic hydrocarbons un- 
der the effect of t he UV radiation field (e.g. lDartois et al.ll2005l; 
iJones et al1ll990h . (ii) the release ("decoagulation") of the VSGs 
and PAHs from the dust aggregates (present in dense clouds) 
at different depths into the PDRs or (iii) a size segregation ef- 
fect due to the grain dynamics induced by the anisot ropic radia- 
tion fie ld, as already mentionned for NGC2023 by Aberg el et alj 
(l2002h . All these processes could be at work with different effi- 
ciencies depending on the depth into the cloud and/or on the size 
and nature of dust. 

This strong evolution of the PAH relative abundance be- 
tween the dense and diffuse medium in PDRs could be a clue for 
the interpretation of the observed (PAH 8 fj.m)/24 //m_and_(PAH 
Sjjm)/ 1 60 jum variations in the SINGS galaxies (e.g. lBendo et alj 
12008 ). Such a PAH relative abundance evolution in the bright IR 
emitters that are PDRs will have an impact on the estimation of 
the star formation activity of galaxies using PAH emission. 
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